// #include "ms_voxelimage.h"
// #include "ms_sphere.h"
// #include "ms_cylinder.h"
// #include "ms_misc.h"

#include "ms_voxelimage.h"
#include "ms_cylinder.h"
#include "ms_misc.h"
#include "unistd.h"
#include "ms_distribution.h"
#include "fe_vecmat.h"

#include "ms_distribution.cxx"
#include "ms_sphere.cxx"
#include "ms_misc.cxx"
#include "ms_voxelimage.cxx"
#include "fe_vecmat.cxx"


#include <stdlib.h>
#include <stdio.h>
#include <iostream>
#include <fstream>
#include <gsl/gsl_histogram.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_statistics.h>

using namespace std;


int main(int argc, char ** argv)
{
	const gsl_rng_type * T;
	gsl_rng * r;
	int i;
	int n = 20;

	double mu = 1.0;
	bool store = true;
	Vector W;
	Vector Blub;

	gsl_rng_env_setup();
	T = gsl_rng_default;
	r = gsl_rng_alloc (T);
	RandomDistribution* useDistribution = new ExponentialDistribution( mu ); 
	for(int j = 0; j<20 ; i++)
	{
		Blub[j] = 0.;
	}
	useDistribution->GetRandomVector( Blub , n );

	Histogram Blob( Blub , W , n , store );

	

	printf ("\n");
	gsl_rng_free (r);
	return 0;
}
